Ageing effects in single particle trajectory averages 
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We study time averages of single particle trajectories in scale free anomalous diffusion processes, 
in which the measurement starts at some time ta > Q after initiation of the process at the time origin, 
t = Q. Using ageing renewal theory we show that for such non-stationary processes a large class of 
observables are affected by a unique ageing function, which is independent of boundary conditions 
or the external forces. We quantify the weakly non-ergodic nature of this process in terms of the 
distribution of time averages and the ergodicity breaking parameter which both explicitly depend on 
the ageing time ta- Consequences for the interpretation of single particle tracking data are discussed. 

PACS numbers: 87.10. Mn,02.50.-r,05.40.-a,05.10.Gg 
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Ergodicity in the Boltzmann sense, the equivalence of 
sufhciently long time and ensemble averages of a physical 
observable, and time translational invariance are hall- 
marks of many classical systems. In contrast, ergodicity 
violation and ageing effects are found in the dynamics 
of various complex systems including glasses [ij , blinking 
quantum dots [l], weakly chaotic maps fsj, and single 
particle tracking experiments in living biological cells [J| . 

As originally pointed out by Bouchaud [ij], when the 
lifetimes t of states of a physical system are power-law 
distributed, ipit) ~ r"/fi+" (0 < a < 1), with an infin- 
ite mean sojourn time (t), both weak ergodicity breaking 
and ageing in the following sense occur: time averages of 
physical observables remain random even in the limit of 
long measurement times and differ from the correspond- 
ing ensemble averages [5-14|. Furthermore ensemble av- 
eraged correlation functions of observables, taken at two 
time instants <2 and ti, are not solely functions of the 
time difference \t2 — ti\. With this breakdown of sta- 
tionarity, the statistical properties of such systems are 
no longer time translation invariant. Consider, for ex- 
ample, a particle undergoing a random walk in a ran- 
dom environment, such that the particle hops a finite 
distance to the left or right with equal probability, but 
sojourn times between jumps are distributed like 'ip{t). 
Such a model, called the continuous time random walk 
model (CTRW), was introduced by Scher and MontroU 
in the context of charge carriers in disordered systems 
[l5| . Its intriguing statistical properties become appar- 
ent when studying time averages, like the time averaged 
mean squared displacement (TAMSD), which is com- 
monly used to analyze the diffusive properties measured 
in single particle tracking assays. Assume that the pro- 
cess starts at time t = Q. From a single trajectory x{t), 
observed in the time interval {ta, ta+T), the TAMSD is 
typically determined through the definition 



with the lag time A and the measurement time T. 

For Brownian motion, in the limit of large T, we ob- 
tain the expected behavior, namely, — 2KiA. In 
this case, the process does not age: the result does not 
depend on the choice of ta- Moreover Brownian mo- 
tion is ergodic: the result for the TAMSD is not ran- 
dom, and the same as found for an ensemble of Brownian 
particles, (a;2(A)) = 2KiA Thus the diffusion con- 
stant Ki can be determined either from an ensemble 
of trajectories as originally done by Perrin 17[ or, as 
conceived by Nordlund 18], from a single particle tra- 
jectory. For anomalous diffusion, defined in terms of 
(x'^iA)) = 2XaA"/r(l -I- a) with the generalized dif- 
fusion constant Ka, the equivalence between time and 
ensemble averages as well as the time translation invari- 
ance generally breaks down. One of our central results is 
that for CTRW processes with long tailed ip{t) the aged 
TAMSD (P) for A < T follows 



Aa^jta/T) ff(A) 

r(l + a) 
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ta 



'x{t + A) -x{t)Y 
T- A 



dt. 
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where Aq(z) = (1 -I- z)" — z", and g is a function of 
A only. Another important finding is that 6^ remains 
a random variable, whose statistical properties explicitly 
depend on both process age ta and measurement time T. 

Previous work revealed deviations from standard er- 
godic statistical mechanics by studying time averages in 
the interval {ta = 0,r) [i-|l2|. Namely, in these works 
the start of the measurement coincides with the start of 
the process jl9j |. However, in experiment, the observed 
particle may be immersed in the medium long before we 
start our observation. In fact in some cases we may not 
even know when the process was initiated. 

In this Letter, we derive the ageing renewal theory for 
CTRW processes and study the dependence of time av- 
erages of physical observables on the starting time ta of a 
measurement after system preparation at t = 0. Accord- 
ing to Eq. the time intervals (0,T) and {ta,ta + T) 
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are not equivalent: In complete contrast to Brownian 
motion, the statistical properties of CTRW trajectories 
will depend on the observation window, the process ex- 
hibits ageing. The remarkable property of Eq. ^ is that 
corrections due to ageing enter in terms of a unique pre- 
factor depending on ta and the measurement time T. We 
call this Aq, the ageing depression function, which is inde- 
pendent of A. This function contains all the information 
on ageing, and is universal since the formula applies for 
any external force or boundary condition, and, as shown 
below, also holds for a large class of physical observables. 
The function g{A) contains the complete A-dependence. 
For instance, we find 5(A) ~ A for free motion ml. Ml and 
g{A) ~ A^~" at long times under confinement [9l. llOj. 

Ageing renewal theory. In a CTRW, the position co- 
ordinate X of a random walker is an accumulation of ran- 
dom jump lengths, x{n) = J2^=o^^i- simplest, 
unbiased version of the model, the Sxi are independent, 
identically distributed (IID) random variables with zero 
mean, and we assume also finite variance a^. Jumps are 
separated by random IID waiting times, drawn from the 
common distribution ■>p{t) T°'/t'^+°', with < a < 1. 
This implicitly defines a counting process n(t), the num- 
ber of steps up to time t. The statistics of the overall 
diffusion process x(t) = x{n(t)) are to be derived from 
the properties of both of its constituents, a method com- 
monly called subordination (20l - |22j . For a large variety 
of physical applications of CTRW and non-aged renewal 
theory see [15|, |23|, ■ We first discuss the counting pro- 
cess n{t)^ emphasizing the role of ageing. 

Since waiting times are independent, the counting pro- 
cess n{t) is a renewal process, which we assume to start 
at t = 0. To study the ageing properties of the system, 
we consider na{ta,t) = n{t+ta) ~~n{ta), the number of re- 
newals in the interval (ta, ta+t). The corresponding prob- 
ability density in double Laplace space, (tail! — > {sa, s), 
in the scaling limit of large times becomes [H, [2^ 



p{na]Sa,s) = 5{na) 



1 



h{Sa,s) 



(3) 



Here the probability of the waiting time for the first jump 
to occur after start of the measurement at ta is 2614281 



h{Sa,s) 



Sa{Sa - S) 



^ h{ta,t) 



t"ita+t) 



(4) 



In the Brownian limit a — 1, the number of jumps n and 
real time t are equivalent, p{na', ta, t) = d{na — t/r). 

This formalism allows for a direct calculation of the 
average of any function of na{ta, t). For instance, the qth 
order moment becomes 129] 



{nl{sa,s)) 



nlp(na;sa,s) dUa 



r(9 + i)--f 



SaiSa - s) Sl 



+aq 



(5) 



After double Laplace inversion, we find 

{nl{ta,t)) = T{q + l)/[T{a)T{l + aq-a)\ 

't + ta 



B 



.1 + aq — a,a 



(6) 



t + ta' 

where B{z;a,b) is the incomplete beta function [sof . 
Thus, the number of steps taken during a time inter- 
val of length t is not stationary in distribution: the 
moments for the period [0,t] are clearly different from 
those for [ta,ta + t]. Indeed we see from Eqs. ^ and 
^ that the process gets slower and eventually stalls as 
ta 00. Concurrently the i-dependence changes with ta'. 
(n?(0,i)) ~ t"i at ta = 0, but {nUta,t)) ~ t^-H^-^+^i 
for large but finite ta/t. 

For the probability density (jS)) we obtain 



p{na;ta,t) 



[t/rr 



H 



ma{t/ta)] 5{na) 

(2 



1,1 



f m„(iAa)r(2 
2a, a) 



(0,1) 



a) 
(7) 



_{t/T)- 

for ta ^ t, in terms of an i?-function 'sil, [s^] and the 
probability to make at least one step during [ta ,ta +t], 

B{[l+ta/t]-^ 



ma{t/ta) = 



,1 



a, a) 



(8) 



r(l - a)T{a) 

Again, we emphasize the explicit dependence on ta. 
When compared to the form for ta = 0, , the most 
striking difference in Eq. ([7]) is the occurrence of the 
(5(na)-term: For a non-aged process {ta = 0), we have 
nia = 1, while an aged process always has a nonzero 
probability of not performing any steps at all within the 
chosen time window, its amplitude approaching 1 algeb- 
raically, m„ - {t/ta)^~°'/[T{a)T{2 - a)], as ta > t. 

Ageing CTRW. We proceed by adding the position co- 
ordinate to our description. For IID jump distances 6xi , 
the jump process x{n) converges in distribution to free 
Brownian motion in the scaling limit of many jumps. The 
anomalous diffusion process x{t) — x{n{t)) inherits the 
ageing properties of the counting process discussed above. 
. To see this, consider the qth order TA moments 



H^;ta,T) 



{\x{ t + A)-xm 

T- A 



dt, (9) 



which are useful to characterize experimental data 
For free Brownian motion x{n) we know that 

2r(<z) 



{\x{n2) - 2^(^1)1'^) 



1^2 - nil 



q/2 



(10) 



r{q/2) 29/2 

Since x{n) and n{t) are independent processes, 
use conditional averaging to evaluate the integrand of 
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Eq. ([9]). From above result 

{\x{t + A)-x{t)\'^) = 

= r(<z + i)/[r(a)r(i + 

X [K^ {t + A)"]"/' B 



IGJ), it follows that 
2r(g) c 



r{q/2) 29/2 
aq/2 - a)] 
A 
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t + A 



,1 - a + aq/2,a j(ll) 
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m =1 



m =0.063 




Figure 1: Time averaged mean squared displacement 
S'^{A;ta,T) for individual free CTRW trajectories (full sym- 
bols) and the averages according to Eq. (|12p (bold black 
lines). Left: Measurement starts at ta ~ 0, so that = 1. 
Right: Aged process, ta — 10^^ (a.u.), so that a large fraction 
1 — TTia ~ 94% of trajectories is suppressed in the log-log plot. 
The parameters are a = 1/2, t — 1, — 1 and T = 10^. 



where we identify = cr^/(2T") The gth order 

TA moment ^ in the limit A ^ T then becomes 



A4ta/T)[i^„A"]'/'r(g+l) /A 



T{a + l)T{2-a + aq/2) \T 



l-Q 



(12) 



which is of the general shape (0) with g(A) /^^-a+aq/'^^ 
Interestingly, we see the special role of the TAMSD [q = 
2), for which the A-scaling is independent of a. 

Fig. [T] shows simulations results for the TAMSD. If 
evaluated during the initial time period, ta — 0, the 
TAMSD for individual trajectories scatter around the en- 
semble average (5^), Eq. ([T2l) . In contrast, for the aged 
process {ta ^ T) the ensemble average (S^) appears much 
lower than the shown individual trajectories. This is due 
to the fact that a significant fraction 1 — nia of particles 
do not move during the entire measurement time. Such 
trajectories are naturally not visible in a logarithmic plot, 
while being included in the calculation of ((5^). 

In a biological cell, the diffusive motion of a tracer 
particle is spatially confined, or a charge carrier in a 
disordered semiconductor experiences a drift force. To 
address such systems we now determine the TAMSD in 
the presence of an external potential. We start with a 
Hookean force —Ax and address the more general case 
below. To that end the continuum approximation of the 
process is represented by the Langevin equation [s^l 



dx/dn = —Xx{n) + ^{n), 



(13) 



where ^{n) is white Gaussian noise with {S,{ni)^{n2)) = 
cr^(5(n2 — rii). In other words, x{n) is a stationary 
Ornstein-Uhlenbeck process. Its increments are Gaus- 
sian variables, characterized by the variance {[x(n2) — 
x{ni)]'^) — a^[l — exp(— A|n2 — ni|)]/A. To calculate the 
TAMSD, we follow the approach outlined for the free 
particle, a general approach to compute this and similar 



quantities will be provided below. The result reads 



S^{A;ta,T)) ^ i^^^^:'^} f,lt Ec.,2 (-A.A") (14) 



r(i 



in terms of the generalized Mittag-Leffler function j3£ 
where A^ = A/r". We deduce the limiting behavior 



(52 



2AMT)K, 
r(l + a)Ti 




r(2-Q) ' 



A « A-^/" 
A » A-^/" 



(15) 



Eq. ([T4|) is another special case of Eq. ([2]) . Interestingly 
the entire dynamics are multiplied by the unique factor 
Aq. Let us now test the generality of this feature. 

Consider the time average of some observable F{x2, xi) 
along the trajectory. 



F{A;ta,T)) = 



{F{x{t + A),xm 
T- A 



dt. (16) 



F may represent moments {F = d'^ , F(x2, xi) = \x2 — 
xi I*) or the TA of a correlation function. We only require 
that the jump process x{n) and the function F fulfill 



{F{x{n2), x{ni))) = f{n2 - ni). 



(17) 



For instance, we have /{n) — a'^n for the second mo- 
ment of unbounded motion [cf. (fTO)) ]. or /(n) ~ cr^[l — 
exp(— An)]/A for the TAMSD in an harmonic potential. 
Condition ([T7)) is fulfilled whenever x{n) is a station- 
ary process (e.g., equilibrated Brownian motion). Al- 
ternatively, one may consider a process with stationary 
increments (e.g., unbounded Brownian motion), when 
F{x2, xi) — F{x2 — xi). In these cases we find 



FiA;ta,T)) = 



ta+T-A 



f{na)p{na]t,A) 



T - A 

where p{na',ta,t) is defined in Eq. ([3]). We obtain 

A4ta/T) 5(A/t) 



duadt, 
(18) 



FiA;ta,T) 



C 



T{l + a) {T/tY- 



(19) 



at short lag times A ^ T, with the constant C = /(O). 
The function g in Laplace space is defined as [1^ 



9{s) 



„2q-2 



^{/(n)-/(0);n^s"}. 



(20) 



Comparing with our specific results (IT^ . and ([T5|), 
we identify all relevant terms in TAs: (i) In the limit 
A — >■ 0, the TA reduces to the constant C, which equals 
the expectation value of the observable when measured 
at identical positions. For example, if we study correla- 
tions in an equilibrated process, F{x2,xi) = X2X1, then 
C ~ (x^) is the thermal value of x^. Conversely, C 
naturally vanishes if we are interested in TA moments 
of displacements, i^(x2,Xi) = |x2 — xij'^, so it did not 
appear previously, (ii) The lag time dependence enters 
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Figure 2: Numerical validation of Eq. (I2ip . for various a, 
boundary conditions, and ta (see key). Each point in the 
graph represents an individual trajectory. Parameters are r — 
1 (a.u.), o-^ = 1, A = 100, T = 2 X lO*^. is either (non- 
aged), or for specific a chosen such that ma ~ 0.054 (aged). 



through the function (/(A). For example, if f{n) ~ n*^, 
then C = 0, and in Laplace space g{s) ^ s"~^~"'', which 
implies 5(A) ~ Ai-"+"'J as in Eq. (iii) The ageing 

depression function Aq only depends on the ratio ta/T 
and the parameter a, and due to a factor T°-~^ any TA 
converges to the constant C as T — ^ 00. Note that this 
dependence on ageing and measurement time ta and T 
is universal in the sense that it is indifferent to the spe- 
cific choice of observable F or model of the jump process 
a;(n), but is directly deduced from the nature of the age- 
ing counting process n(t). Also note that in the Brownian 
limit a = 1, Eq. dUl) reduces to {F{A;ta,T)) = /(A/r), 
restoring the ergodic equivalence of ensemble and time 
averages, and the stationarity of the process. 

Distribution of TAMSD. Due to the scale-free nature of 
the distribution tp{t) of waiting times all TAs of physical 
observables, e.g. S'^, remain random quantities, however, 
with a limiting distribution for the dimensionless 
ratio ^ = (52/((5^) d, [13, Hlj. As contributions to time 
averages of the form ^ occur at time instants when the 
particle performs a jump, we expect that in the sense of 
distributions both 6^ and Ua should be equivalent, = 
cUa, for some non-random, positive c. In other words. 



S^A;ta,T) d na{ta,T) 



{S^A;ta,T)) {na{ta,T))' 



(21) 



for A ^ T. We may thus deduce the statistics directly 
from the underlying counting process. Fig. [2] provides 
numerical evidence for this argument in the case of a free 
particle and a particle in a box for several values of a. 

The distribution 0(f) for = is related to a one- 
sided Levy stable law [6|. In the opposite case ta ^ T, 
combination of Eqs. ((2T|) and ([7]), yields 



m ~ [l~m^{T/ta)]5{^) 

l-a 



{T/ta 



r(a) 



-H. 



1,0 
1,1 



.{T/ta^ 



m^{T/ta)T{2 - 

(2 -2a, a) 
(0,1) 



a) 



.(22) 




0.21 



0.04 



0.02 





— Eq. (22) 
. a = 1/2 

a = 3/4 







Figure 3: Scatter density <^(^) for different ct and m,a, see 
text. Lines: Eq. (7) from ;6i (Left) and Eq. (I22j (Right). 
Symbols: Simulations of free CTRW. Note that the area under 
the curves for the aged process in the right panel is not unity, 
since the fraction 1 — rria of immobile events is not shown. 
Same parameters as in Fig. [S] 
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Figure 4: Ergodicity breaking parameter (|23|) as function of 
a (Left) and ta/T (Right). Notice that the non-ergodic fluc- 
tuations become larger with increasing ta- 



for A < T < ta- In 0(f), ma{T/ta) is the weight of the 
continuous part. The probability for not moving dur- 
ing the whole measurement period, (f = 0), approaches 
unity as ~ {T /taY~°' ■ If conditioned to measurements 
with f > 0, we also find the scaling f ^ {ta/Ty~°'. In 
Fig. |3] we demonstrate excellent agreement of Eq. (|22|) 
with numerical simulations. 

Deviations from ergodic behavior are quantified by the 
ergodicity breaking parameter 6] , for which we obtain 



EB 



<52 



B{[l + ta/T]-\l + a,a) 
[I -{1 + T/ta)-"]" 



(23) 

depending only on the ratio ta/T. At = EB re- 
duces to the result of Ref. [6j, while for ta ^ T, we 
find EB - 2(ta/T)i-"/[a(l + a)]. In the non-aged case 
EB is bounded, < EB < 1, depending on the value 
of a only. In contrast we find that EB may diverge in 
the limit ta/T 00. This implies that the non ergodic 
fluctuations are much larger in the aged regime under 
investigation. We show the behavior of EB in Fig. ID 

Conclusions. We investigated the effects of ageing on 
TAs of physical observables. Previous calculations of TAs 
tacitly neglect the fact that often the preparation of the 
system and start of the measurement do not coincide. 
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While this does not cause any problems for ergodic sys- 
tems with rapid memory loss of the initial conditions, 
in general this cannot be taken for granted in processes 
of anomalous diffusion. Here we showed for the case of 
CTRW with diverging characteristic waiting time that 
TAs of arbitrary physical observables carry the common 
factor Aq. This ageing depression function is universal 
in the sense that it only depends on the process age ta 
and the measurement time T. All details such as con- 
finement effects enter through a single function, 5(A). 
The structure of this result was shown to hold for a large 
class of physical observables. We also see that the ageing 
of the process has a pronounced statistical effect, split- 
ting the population into two: the mobile fraction TOq 
and the immobile one whose amplitude 1 — ma grows 
with ta/T. Knowledge of this effect is significant for the 
quantitative physical interpretation of experimental data. 
Finally, since renewal theory is applicable to many sys- 
tems, our results with minor changes should be relevant 
more generally, for instance, to the Aaronson-Darling- 
Kac theorem in infinite ergodic theory or for counting 
the number of renewals in blinking quantum dots. 

We acknowledge funding from the Compint graduate 
school, the Academy of Finland (FiDiPro scheme), and 
the Israel Science Foundation. 
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